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Doin' the twist: Secular changes in the surface differential 
rotation on AB Doradus 
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ABSTRACT 

We present measurements of the rotation rates of individual starspots on the rapidly 
rotating young KO dwarf AB Doradus, at six epochs between 1988 December and 1996 
December. The equatorial rotation period of the star decreased from 0.5137 to 0.5129 
days between 1988 December and 1992 January. It then increased steadily, attaining a 
value of 0.5133 days by 1996 December. The latitude dependence of the rotation rate 
mirrored the changes in the equatorial rotation rate. The beat period between the 
equatorial and polar rotation periods dropped from 140 days to 70 days initially, then 
rose steadily. The most rigid rotation, in 1988 December, occurred when the starspot 
coverage was at a maximum. The time-dependent part of the differential rotation is 
found to have Af2/f2 ~ 0.004, which should alter the oblateness of the star enough 
to explain the per iod changes observed in several close binaries via the mechanism of 
Applegate (1992)[ 
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1 INTRODUCTION 

The back-reaction of Lorentz forces on the fluid circulation 



rotating K0 dwarf AB Doradus, providing independent sup- 
port for this theory. We apply a new method for determi ning 



in the cpiivccLivt! zone of a luLaliiig laLe-Lype sLai is uf fuiida- 
mental Umpuitaiice in theories of sl i dl 
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i magnetic-field gun- 

eration. Over the last decade or so, a consensus has emerged 
among dynamo theorists that this back-reaction could have 
a significant influence on the differential rotation pattern 
in the co nvective zones of the Sun and active late-type stars 



the latitudes and rotation rates of ind ividual starspots (Col- 



lier Cameron, Donati & Semel 2002) to six sets of archival 



(see, e.g. Brandenburg et al 1991; Moss et al. 1995; Kitchati- 



nov et al. 199S ) 



The quasi-cyclic orbital period changes observed among 
short-period binaries with magnetically-active components 
provide further indirect evidence that substantial changes in 
differential rotation co uld be taking pla c e on active binary 
components. Work by Aj)pl£gate^l^92j , [La nza, Rodono & 



Rosner (1998)| and Lanza fc Rodono (1999) ' suggests that 



small, magnetically-modulated changes in stellar differential 
rotation will alter the gravitational quadrupole moment of 
the active star. This can be sufficient to produce orbital 
period changes as large as the fractiona l 8P/P ~ 1.5 x 10 -4 



time-resolved echelle spectroscopy of AB Dor, spanning the 
period from 1988 December to 1996 December. 



2 OBSERVATIONS AND DATA REDUCTION 

The details of the instruments and observing procedures 
used to secure the six data sets have been published else- 
where, so we give only a brief summary in Table [j]. 

The observations from the three earliest runs were re- 
extracted from the raw data using the Starlink ECHOMOP 
optimal extraction routines, to ensure consistency with the 
later data. The spectra from all 6 years included the Ha re- 
gion, in which numerous weak, narrow telluric lines of H2O 
and O2 are present. We used spectral subtraction (Collier 



Cameron et al. 2001) to isolate the travelling distortions pro- 



observed in HR 1099 by ponati (1999) . 

In this paper we report observations of long-term 
changes in surface differential rotation on the young, rapidly 
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duced by sta rspots in the mean profile, and least-squares de- 
convolution (Donati et al. 1997) to stack up the residual pro- 



file information in the large number of known photospheric 
lines recorded in the echellograms. The number N of lines 
used in each year ranges from 140 to 2000, and is given in 
Table [l], column 3. Since the spectra in all years were ex- 
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Table 1. Summary of observing runs. Details of the observations are given in the papers listed. The phase ranges in the fourth column 
are those in which the complete signatures of individual spots could be observed on at least two nights of the run observing concerned. 

References 



UT Dates 


Telescope /Instrument 


No of 


Phase coverage 






lines 




1988 Dec 16, 19 


3.9-m AAT/UCLES 


139 


0.09:0.40 


1988 Dec 21, 23 


3.6-m ESO/CASPEC 


326 


0.0:0.38, 0.94:1.00 


1992 Jan 18, 19, 20 


3.9-m AAT/UCLES 


1134 


0.08:0.62 


1993 Nov 23, 24, 25 


3.9-m AAT/UCLES 


1566 


0.40:0.85 


1994 Nov 15, 16, 17 


4.0-m CTIO/echelle 


619 


0.04:0.18, 0.28:0.63 


1995 Dec 07, 11 


3.9-m AAT/UCLES 


1936 


0.55:0.95 


1996 Dec 23, 25, 27, 29 


3.9-m AAT/UCLES 


1964 


0.50:1.00 



3ollier Cameron ct al. (1990) 




Dollicr Cameron ct al. (1990' 




ZJollicr Cameron & Uuruh (1991) 




Umuh. Collier Cameron & Cutispoto (1995) 


Collier Cameron ct al. (1999) 




Donati & Collier Cameron (1997) 




Donati ct al. (1999) 





posed to a signal-to-noise (S:N) ratio of 100 to 120, the S:N 
of the deconvolved profiles scales approximately with y/N. 

The deconvolved profiles were placed in the heliocentric 
reference frame, to an accuracy better than 100 m s _1 , using 
the telluric lines as velocity references during the deconvolu- 
tion procedure. The resulting time-series of residual profiles 
were subjected to a matched-filter analysis, yielding mea- 
sures of the spot area, radial velocity amplitude about the 
stellar centre of mass, and rotation period, together with es- 
timates of their uncertainties. More detailed descriptions of 
the deconvolution, spectral subtra ction and matched-filter 



analysis procedures are given by Collier Cameron et al 



(2001)Qand |Collier Cameron, Donati fc Scmel (2002) 



As 

in this earlier paper, we used a limb-darkening coefficient 
u — 0.77 and vsiw, 
matched filter. 



91 km s in constructing the 



3 RESULTS 

The rotation periods of the spots detected in each season's 
data are plotted against their rotational velocity amplitudes 
K — Q(9)R+ cos f? sin i in the left-hand panels of Figure |l| 
Here Q(0) is the rotation rate at latitude 9, and J?* is the 
stellar radius. 

In the 1995 and 1996 seasons, outlying points caused 
by aliasing between closely-spa ced pairs of spots observed 
several day s apart can be seen (Collier Cameron, Donati & 
Semel 2002). These were excluded from the next stage of the 
analysis. 

The stellar radial velocity v r ~ 32.5 km s _1 and 
vsini = 91 km s _1 are both uncertain by ±1 km s _1 or 
so. The systematic errors in the derived quantities K and 
P due to these uncertainties are smaller than the la error 
bars sho wn in Figure | l| 



As Collier Cameron, Donati & Semel (2002) noted, 



there is an intrinsic scatter of the individual spot rotation 
rates about the mean differential rotation pattern. The rms 
magnitude of this scatter is equivalent to an additional un- 
certainty of ±0.00017 day in the rotation period of each spot. 
After adding this error in quadrature to the measured uncer- 
tainty in each spot's period, we fitted a differential rotation 
law of the form 



f2(0) = ^equator ~ n beat sin2 



(1) 



to each season's data. The optimal solutions for the equa- 
torial rotation period ^equator = ^ /^eqa&tor and the 
equator-pole beat period Ajgat = ^/^beat are suown to- 
gether with their 68.3% and 95.4% confidence regions in the 



Table 2. Differential rotation parameters derived from 2- 
parameter fits to the rotation rates of individual spots in 1988 
December, 1992 January, 1993 November, 1994 November, 1995 
December and 1996 December. The second column gives the num- 
ber of spots contributing to each fit. 



ID 


No of 
spots 


^equator 
(days) 


^"beat 
(days) 


x 2 


47517.1 


12 


0.51367 ± 0.00010 


13611? 


32.4 


48640.9 


12 


0.51285 ±0.00012 


69+1° 


11.2 


49316.3 


15 


0.51290 ±0.00010 


71+ 6 

' -6 


47.1 


49672.5 


13 


0.51299 ± 0.00020 


94l 3 2 6 


55.2 


50061.2 


16 


0.51292 ± 0.00008 


88+^ 


15.9 


50443.7 


19 


0.51326 ±0.00007 


109+9 


53.3 



right-hand panels of Fig. The inner contour is the locus 

X 2 = Xmin + whose extremities in P e quator and p beat 
give the one-dimensional la error bars on the two periods. 
The outer contours, at \ 2 ~ X^in ± 2.3 and 6.17, contain 
68.3% and 95.4% of the joint probability respectively. 

The fits to the six years' data are summarized in Table ^ 
and Fig. ^. The equator-pole beat period of AB Dor displays 
significant variability from one observing season to the next. 
Between 1988 December and 1992 January the differential 
rotation rate in the star's surface layers doubled, with the 
beat period falling from 140 to 70 days. The equatorial rate 
increased and the high-latitude rotation rate dropped. This 
increased shear pattern persisted into 1993. In 1994 and 1995 
the rotation rate at latitude 40° and northward increased 
while the equator continued to rotate rapidly. In 1996 the 
high latitudes continued spinning up as the equator began 
to show signs of slowing down. The smallest variation in 
rotation period is seen at K ~ 70 km s _1 (latitude ~ 40°), 
but nowhere on the star is the rotation rate constant. 



4 DISCUSSION AND CONCLUSIONS 

We have found that the differential rotation rate of the 
young K0 dwarf AB Dor varies by a factor 2 on a timescale 
of a few years. While the 8-year timespan of the observa- 
tions presented here is too short to determine whether we 
are seeing a cyclic phenomenon, the apparent smoothness of 
the increase and decline in the differential rotation suggests 
that the changes occur on a timescale of order a decade or 
two. The variation resembles a wave of excess rotation, orig- 
inating at the equator in 1 992 and propagating poleward 



Applegate (1992) and Lanza, Rodono & Rosner (1998) 



have shown that changes in magnetic flux threading the con- 
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Figure 1. The left-hand panels show rotation period P versus rotational velocity amplitude K for candidate starspots in each year 
of observation. "Threshold" denotes they 2 threshold below which spurious spots have been rejected. "Phi" denotes the phase range 
searched, as given in column 4 of Table hi The fitted differential rotation curve is shown for each year. The right-hand panels show the 
contours A\ 2 = 1.0, 2.3 and 6.17 as a function of equatorial rotation period and equator-pole lap time, both in days. 
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Dor 1994 threshold = 25 phi = 0.04-0. ' 8 and 0.28-0.63 
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Figure 1 - continued 
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AO. = 0.048 radian d 1 in the equator-pole beat frequency, 
we find that AQ./Q, ~ 0.004. 

Such a large modulation of the surface differential ro- 
tation will alter the star's oblateness sufficiently that, if AB 
Dor were in a close binary system, it would be expected 
to produce observable long-term orbital period changes. By 
way of comparison, the K2V primary of the binary system 
V471 Tau has a spectr al type and rotatio n period very simi- 
lar to those of AB Dor. Applegate (f992) points out that the 

" modulation of the orbital 



observed 20-year, A P/P ~ 10 
period in V 471 Tau ( fckillman 



47500 48000 48500 49000 49500 50000 50500 — aL 1994 l) would require 



Patterson 1988; [banoglu 



Julian Date-2400000.0 

(a) 




47500 48000 48500 49000 49500 50000 50500 
Julian Date-2400000.0 

(b) 

Figure 2. The equator-pole beat period Pfo ea t (upper panel) and 
equatorial rotation period Peq (lower panel) versus Julian date. 
The data points are shown for (left to right) 1988 December, 
1992 January, 1993 November, 1994 November, 1995 December 
and 1996 December. 



vective zone during a stellar magnetic cycle may alter the 
viscous transport of angular momentum sufficiently to pro- 
duce substantial changes in differential rotation. AB Dor 
was rotating most rigidly in late 1988, at a time when the 
star's maximum light level was depressed by 0.2 mag or so 
relativ e to photometry se cured in the early 1980s and mid- 
1990s (Kiirster et al. 1997). The increase in differential rota- 
tion between 1988 and early 1992 coincides with a 0.1-mag 
rise in the maximum brightness of the star. The maximum 
brightness of the star remained steady between 1992 and late 
1996, as the differential rotation began gradually to decrease 
again. The timescales of the long-term changes in starspot 
coverage and differential rotation thus appear to be compa- 
rable, though we cannot yet determine how closely the spot 
coverage and differential rotation rate are linked. 

These departures from the time-averaged differential ro- 
tation pattern are between one and two orders of magnitude 
greater than the waves of excess and deficit rotational veloc- 
ity - sometimes described as "torsional oscillations" - that 



migrate equatorward during the solar spot cycle (Howard 



Labonte 1980; Labonte & Howard 1982). These depar- 



tures from the mean solar rotational velocity have ampli- 
tudes of order 5 m s _1 . At the equator of AB Dor, however, 
a change in period from 0.5137 to 0.5128 days corresponds to 
a 200 m s _1 change in vsini. If we consider the full change 



variable differential rotation in 
the K star, with Afi/O = 0.0032. The observed changes in 
AB Dor's differentia l rotation are of precisely the magni- 
tude required by the Applegate (1992) model to explain the 
orbital period variations in V471 Tau. 
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